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IV 


1.  Introduction  and  Definition 


This  report  is  a  further  development  of  previous  reports  (1-4)  based  on  the  concept  of  the 
operational  equations  of  state.  It  is  suggested  that  the  APA-  equation  of  state  (EOS)  depends  on 
three  functions  of  one  variable  each  that  can  be  presented  in  the  following  simple  form: 

S(V,E)  =  A(V)+  f  ,  (!) 

'  ’  K  '  l  P(4)+Mv) 

Thus,  in  accordance  with  the  general  paradigm  of  the  operational  EOS,  the  APA -EOS  is  the 
integral  operator  transforming  the  functions  A(V),  P(£) ,  and  A (V) into  the  third  function  of 

two  variables,  called  entropy  S  (V,  £) . 

In  their  own  turn,  the  functions  A(V),  P(E') ,  and  A(V)  are  also  certain  operators  of 

experimental  data.  Eventually,  the  APA -EOS  can  be  treated  as  the  operators  transforming 
experimental  measurements  into  EOS  rather  than  algebraic  functions. 

The  postulated  EOS  equation  1  is  fixed  in  terms  of  the  functions  A (V) ,  P(£) ,  and  A(V) .  We 

call  it  the  skeleton  of  the  operational  EOS.  However,  our  goal  is  to  present  the  operational  EOS 
as  an  operator  whose  arguments  are  not  abstract  functions  but  experimentally  measured  data. 
Those  ultimate  operational  EOS  depend  upon  the  measuring  tools  available.  By  using  the  same 
skeleton  but  different  measuring  tools  and  procedures  the  ultimate  operational  EOS  will  be 
different.  In  particular,  the  ultimate  APA -EOS  depends  on  the  specific  types  of  the  available 
measurements-isochoric,  adiabatic,  isothermal,  etc.  The  general  intermediate  structure 
equation  1  guarantees  that  with  any  options  of  the  chosen  experiments  and  measuring  devices, 
the  ultimate  APA -EOS  will  be  complete  and  thermodynamically  consistent. 

We  remind  the  reader  that  completeness  means  that  the  EOS  contains  the  full  thermodynamic 
information  about  the  substance  and  any  required  thermodynamic  potential  or  function  can  be 
extracted  by  standard  algebraic  and  calculus  operations  and  without  any  additional  physical 
measurements.  Thermodynamic  consistency  means  that  the  APA— EOS  is  totally  consistent  with 
the  first  and  second  laws  of  thermodynamics.  Remember  that  the  most  popular  and  widespread 
EOS  in  the  high-pressure  and  shock-wave  analysis  community,  the  Mie-Gruneisen  EOS,  is 
neither  complete  nor  necessarily  thermodynamically  consistent  with  other  assumptions. 

Generally  speaking,  the  APA— EOS  is  incompatible  with  the  Mie-Gruneizen  EOS  (1,  2). 
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Roughly  speaking,  there  are  no  constraints  in  the  choice  of  the  function  S  =  S  (V,  E)  from  the 

standpoint  of  thermodynamic  consistency.  In  other  words,  there  are  no  modeling  assumptions 
when  working  with  the  most  general  function  S  =  S  (V,  E) .  At  the  same  time,  the  skeleton 

equation  1  imposes  certain  constraints  on  the  substances  under  consideration.  The  assumptions 
can  be  dictated  by  mathematical  or  computational  conveniences  or  generalization  of  a  simple 
classical  model.  The  assumptions  can  also  be  triggered  by  theoretical  or  experimental  physical 
analysis.  In  all  these  cases  of  modeling  assumptions  there  is  the  risk  of  thermodynamic 
inconsistency,  insufficiency,  or  redundancy  of  generated  models.  The  risk  is  greater  when 
dealing  with  the  sophisticated  mathematical  forms  like  equation  1 . 

The  assumption  equation  1  has  a  clear  physical  formulation  as  it  implies  that  the  specific  heat 
capacity  C  at  constant  volume  V  is  a  function  of  the  internal  energy  density  E  only: 

C=C(E ) .  The  opposite  statement  is  also  valid:  assuming  that  the  heat  capacity  is  a  function  of 

the  internal  energy  only,  one  can  deduce  that  the  specific  entropy  S  can  be  presented  in  the 
general  form  equation  1.  From  that  standpoint,  the  APA-EOS  reminds  the  classical  Boyle- 
Mariotte-Gay-Lussac  model  (when  the  heat  capacity  must  be  function  of  the  absolute 
temperature  T  only)  and  the  Mie-Gruneisen  EOS  (when  the  heat  capacity  must  be  a  function  of 
the  entropy  S  only).  At  the  same  time,  the  APA-EOS  obviously  differs  from  both  previously 
mentioned  classical  models  with  one  important  exception.  This  exception  concerns  the  model 
with  the  constant  heat  capacity  coined  in  ( 4 )  as  the  Dulong-Petit  model.  In  this  case,  all  four 
models  are  equivalent  to  each  other. 


2.  Thermodynamic  Identities  for  the  APA  EOS 


By  applying  standard  thermodynamic  identities  to  the  definition  equation  1  of  the  APA-EOS, 
we  get  the  following  formulas  for  the  absolute  temperature  T ,  pressure  P ,  and  heat  capacity  Cv : 


T(V,£)  =  P(£)  +  A(V)  . 


(2) 


P(V,E)  =  [P(£)  +  A(y)] 


A'(V)-A  '{V)]dt; 


l  [P(£)  +  A(V)]  J 


c^v^  =  V{e)  =  Cv{e)- 


(3) 

(4) 


According  to  equation  2,  the  absolute  temperature  for  the  APA-EOS  appears  to  be  the  sum  of 
the  function  P(£)  of  the  internal  energy  density  E  and  of  the  function  A(V)of  the  specific 

volume  V  .  Equation  4  also  states  that  the  specific  heat  capacity  function  Cv  appears  to  be  the 
function  of  the  specific  internal  energy  E  only. 
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(5) 


The  Gruneisen  function  V{V ,E)  for  the  APA-EOS  is  given  by  the  following  formula: 


r(v,E) 

V 


P'(E)A'(V) 


1  A '(v)P'(g)  ,  rd.  W) 

p'(£-)p(£’)+A(V)  J  p(f)  +  A'(V)[p'(£)]2  ' 


Combining  equations  4  and  5,  we  get 

r(V,E)  A'(V)  Cy(E')  A'(V) _ 1  f  j'Cr(f)A’(V) 

V  CV(E)  CV(E)  ?(£•)+ A(V)  CV(E){  *p(#)+A(V)  ' 


(6) 


3.  Recovery  of  the  A,P,  A  Functions  From  Physical  Experiments 


3.1  Recovery  of  the  Function  P(Zs')  From  Isochoric  Heating 

The  way  of  recovery  of  the  function  P(F)  is  quite  straightforward.  To  that  end,  the  isochoric 

heating  is  the  most  relevant  source  of  information.  In  the  following,  all  the  parameters  relating  to 
the  ambient  state  will  be  supplied  with  the  asterisk  mark. 

Let  T*  =  Ty  (  E)  be  the  temperature  dependence  upon  the  internal  energy  density  at  the  ambient 

volume  V  =  V* .  By  definition,  the  internal  energy  change  in  the  isochoric  process  is  equal  to  the 
heat  supply  Q .  Thus,  the  determination  of  the  function  T*  =T*(E )  requires  a  thermometer  and 

a  device  measuring  the  heat  supply.  Then,  equation  2  allows  us  to  recover  the  function  P(E') : 

P  (E)  =  Tv*  (F)  -  A  ( V * )  =  Tv*  (E)  -  A* .  (7) 


We  can  now  rewrite  the  general  APA-EOS  1  as  follows: 


E 

s(y,£)  =  A(y)+J 

E, 


cU; 

%  (£)+A(y)-A*  ’ 


(8) 


and  equations  2  and  3  as  follows: 

t(v,e)=t;(e)+ A(y)-A*  (9) 


and 


p(y,£)  =  [r;(£)+A(y)-A*] 


a 

A'(y)-A'J 


[Tv  (40+A(y)-A*] 


(10) 
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3.2  Recovery  of  the  Functions  A(V )  and  A(V)  From  the  Quasi-Static  Ambient  Isotherm 
3.2.1  Formula  of  the  Potential 

Let  E*  =E*(V )  be  the  equation  of  the  ambient  isotherm  T  =  T  .  Then,  equation  9  gives  the 
following  formula  of  A(V) : 

A(y)-A*=r-r;(4(y)).  (li) 


Substituting  equation  1 1  in  equation  8,  we  get 

E 

S(V,E)  =  A(F)+J 

E, 


_ dl _ 

t;^)+T'-T;(e‘t(V)Y 


(12) 


Let  Sj  =  Sr  (V)  be  the  values  of  the  entropy  on  the  ambient  isotherm  E*T  =  E*T  (V  ) .  Then  we  can 
rewrite  equation  12  as  follows: 


A(v)=s;(v). 


J 


cU; 


T;{t)+r 


(13) 


Inserting  equation  13  in  equation  12,  we  get 


a 

S(V,E)  =  S*T(V)+  J 


Et(V) 


r;{z)+r -t;(e;(v))' 


This  formula  should  be  augmented  with 


(14) 


E*  (V )  -  E*  =  -  J  darP*  (m)  +  T*  [s;  (V )  -  S*  ] , 

V. 

implied  by 

dEdv)  P‘(v)+r<E) 

dV  ry  '  dV 

3.2.2  Temperature  and  Pressure  Formulas 

Using  equation  14,  we  get  the  following  relationships  for  the  absolute  temperature  and  pressure: 

t(v,e)=t;(e)+t* -t;(e;(v)).  as) 
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P(V,E)  _dS*(V)  1  dE*(V ) 

Tv*(E)  +  T*-Tv*(e*(V))~  dV  T*  dV  + 

f  Je  dT;(E‘T(V))/dV 


I  dt 


3.3  Recovery  of  the  Functions  A(V)  and  A(V)  From  the  Quasi-Static  Ambient  Adiabata 
3.3.1  Formula  of  the  Potential 

Let  E*s  =E*S(V)  be  the  equation  of  the  ambient  adiabata  S  =  S*  in  the  (V,E)  space.  Then,  we 
get 


S(V,E)  =  S‘+  \ 


r;(#)  +  A(V)-A" 


Let  T*  =  T*  (V )  be  the  values  of  temperature,  corresponding  to  the  ambient  adiabata  S  =  S* . 
Then,  using  equation  9,  we  get 

i(v)-x =t;(v)-t;(e:s(v)) 

Inserting  equation  18  into  equation  17,  we  get 


X(V,£)  =  S'+  J 


3.3.2  Temperature  and  Pressure 

Using  equation  19,  we  get 

t(v,e) = t;(e)+t ■;  (v)-t;  (e;  (v))  , 

and 


J  d{ 


t;(v)+t;(e)-t;(eh 

yj)dE's(V ) 

t;(v) 

t;(e)+t;(v)-t;(e;(v) 

j  d 

dV 

[r;(F)-7-;(£;(F))] 

E(z)+t;(v)-t;(e;(v)) 

T 

dV 
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3.4  Recovery  of  the  Functions  A(V )  and  A(V)  From  the  Quasi-Static  Ambient  Isotherm 
and  the  Hugoniot  Adiabata 

3.4.1  Formula  of  the  Potential 

When  dealing  with  this  sort  of  available  measurements  we  can  recover  the  function  /l  ( V )  as 
noted  before  and  deal  further  with  the  equation  containing  the  unknown  function  oc(  V)  only: 

E 

S(V,£)  =  A(V)+  \deK{V,e),  (22) 

E, 


where 


k(v,e)  = 


T;(E)+r-T;(E;(v)y 

It  is  important  to  realize  that  K (V,  £)  should  be  treated  as  the  already  recovered  function. 


(23) 


Using  equation  22,  we  get  the  following  formulas  of  the  absolute  temperature  and  pressure: 

wv) _ 1  hceK{V’C} 

K(V,E)  (  >  K(V,E)l 


t(v’e)=it7 hr,-p(v.E)=-  1 


SV 


(24) 


K(V,E ) 

Let  EH(V)  be  the  energy  measurements  on  the  Hugoniot  adiabata.  Then,  using  the  equation  of 
the  Hugoniot  adiabata  we  arrive  at  the  following  formula  for  the  function  A(V) : 

,  .  ,  f  ,  ,  ^2EH(n)-2E.  +  P.(n-V.)  vr  M7)  dK(n,s) 

A(V)  =  A  +jd7JK(T?,EH(r?))  aK,) -  W - ,—  \dr1  j  de  ^  ,  (25) 

TJ  y *  ^ 


drj 


which  allows  us  the  present  the  operational  EOS  for  hydrocode  in  the  following  form: 


n, 

S(V,E)  =  S*+  J  dsK(V,s)  + 


MO 


v 

jdrjK(ri,EH(fj)) 


2Eh  (; 7 ) - 2 E.  +  P.  (»y-V.)  +  dEH(rj) 


TJ-V, . 


drj 


(26) 


3.4.2  Temperature  and  Pressure 

The  formula  in  equation  26  implies  two  formulas: 

1 


T{y'E)  =  K(^E)  =  T;  (E)  +  r  (V))' 


(27) 
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and 


P(V,E)=  |  de  /  (28) 

Jrv)  S1/  K(v-E) 


V-V. 


3.5  Recovery  of  the  Functions  A(V)  and  A(V)  From  the  Hugoniot  Energy  and 
Temperature  Measurements 

3.5.1  Formula  of  the  Potential 

Let  Eh  (V )  and  Tn  (V )  be  the  energy  and  temperature  measurements  on  the  Hugoniot  adiabata. 
Then,  we  arrive  at  the  following  EOS: 

E 


S(V,E)  =  A(V)+\dsKH(V,s),KH(V,E)  =  - 


1 


t;(e)+th(v)-t;(eh(v)) 


(29) 


We  arrive  at  the  situation  considered  previously.  Therefore,  the  formula  in  equation  26  should  be 
modified  as  follows: 


or  else 


S(V,E)  =  S*+  |  dsK„(V,e)  + 

Eh(V) 


V 

\dTjKH(Tj,EH(rj )) 


2 Eh  (77)  -  2Et  +  Pt  (77  - V, )  dEH  (77) 


TJ-V* 


dtj 


(30) 


S(F,£)  =  S*+  f  ds-rrr - - - 

E.\n  Tv(s)  +  Th(V)-Tv(E„(V)) 


+ 


2 Eh  (77)  -  2Et  +  Pt  (77  -  V* )  dEH  (77) 


rj-V, 


drj 


(31) 


3.5.2  Temperature  and  Pressure 

For  temperature  and  pressure,  we  have  the  following  recovery  formulas: 

1 


T(V,E)  = 


kh(v,e ) 


=t;(e)+th(v)-t;(eh(v)). 


(32) 
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and 


•(V,E)=  ) 


ds 


dKa(V,e)  |  KH{y,EH(V))2EH(y)-2E,+P.(V-V.) 


Eh(V) 


Kh(V,E)  8V 


kh(v,e) 


v-v. 


(33) 


respectively. 


4.  The  Constant  Heat  Capacity  Case 


This  case  was  considered  earlier  in  (4),  using  a  different  approach.  Here  this  important  case  is 
analyzed  from  the  more  general  point  of  view. 

4.1  One  Important  Relationship 

We  begin  with  the  relationship 

f 


p(4)+Hv) 

In  \p(E)  +  X(V)\  ln[^(R)  +  A(V')1 


P’(E) 

implied  by  the  following  chain: 


p(E) 


\n\pU)  +  X{V)~\  In [/?(£’»)  + /I (y)l  f  p"(£)  r  .  ,  .  xl 

m  -  V(E,)  ■ 

We  can  rewrite  equation  34  as  follows: 

J  M.  =ln  +  +J^4^1n[^)+A(V)]. 

E,  p{$)+My)  [p(^)+/i(y)]y(£*)  &  p 

When  p'(E )  =  ME  +  N ,  where  M  and  N  are  constants,  using  equation  34  we  get 


(34) 


(35) 
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(36) 


1 


pW)+My) 


1  ,n  e+a(v) 

M  E.+A(V)’ 


as  it  follows  from  the  chain 

|  df  1  lnME  +  N+MV)  1  [n  E  +  i(V) 
l  p(4)+X(V)  M  ME,  +  N+X(V)  M  E.+A(V)' 


By  inserting  equation  36  in  equation  1,  we  get  the  particularly  simple  formula  of  the  complete 
EOS: 


S(V,E) 


=  a(V)  +  —  In 
v  ’  M 


E  +  X(V) 

E.+X(y)' 


(37) 


or  else 


X(V,£)  =  ff(F)  +  iln[£  +  r(V)]. 


(38) 


With  the  help  of  the  thermodynamic  identities  for  S(V,E ) ,  the  relationship  (equation  38)  allows 
us  to  get  the  relationships  for  the  absolute  temperature  and  pressure 


and 


t(v,e)=m[e+X(v)], 

(39) 

P(V,E)  =  Ma'(V)E+Ma'(V)X(V)+X'(V), 

(40) 

the  last  of  which  is  implied  by  the  chain: 


p(v,e)=m[e+X(v) 


a'(V)  +  - 


i'(v) 


m\_e+X(v)]  J 

M[E  +  X{V)\d'{V)  +  X,(V)  =  Ma'(V)E  +  Md'{y)X{V)  +  X'(V). 
In  the  case  of  constant  heat  capacity,  equations  2  and  4  imply 

n(E)  =  ^(E)=^driC^  =  -(E-E.)  =  T*(E)-X*  , 

C 


(41) 


(42) 


whereas  equation  1  reads 


E 

S(V,£)  =  a(V)  +  cJ 

E-, 


d % 


£-Et.+CJl(V) 


E-Et  +  CA(V ) 

=  a(V)  +  C In - -  -  v  -  . 


CA(V) 


The  thermodynamic  identities  1  and  3  can  be  reduced  to  the  form 


(43) 
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(44) 


T(V,E)  =  ±(E-E,)  +  A.(V)  =  -  E‘+^(V'1  =tl(v,e), 
P(V,E)  =  a’(V)TL(V,E)-^(E-E,), 


W) 


(45) 


the  last  of  which  can  be  also  rewritten  as  follows: 

P(V,E)=(E-Et) 


c  v  ’  x(v) 


-a'(V)A(V). 


(46) 


4.2  Recovery  of  the  Functions  a(V)  and  A,(V)  From  the  Quasi-Static  Ambient  Isotherm 

This  case  in  general  is  described  by  the  relationship  14.  Using  equation  42,  we  can  rewrite 
equation  15  as  follows: 

E-E*(V)  +  CT* 


S(V,£)  =  Rein¬ 


er 


(47) 


Equation  43  implies  the  following  recovery  equations  for  the  absolute  temperature, 

E-E*(V)  +  CT* 


t(v,e)  = 


c 


(48) 


and  pressure 


E-EZ(v)+cr  dS'T(V)  dE^jv) 
'  ’  ’  C  dV  dV 


(49) 


4.3  Recovery  of  the  Functions  a(V)  and  A,(V)  From  the  Quasi-Static  Ambient  Adiabat 

This  data  is  recovered  by  the  formula  equation  19.  Combining  it  with  equation  42,  we  get  the 
following  complete  EOS: 

e-e‘s(v)+ct;(v) 


S(V,£')  =  S'+Cln- 


crs(v) 


(50) 


The  formula  (equation  50)  implies  the  following  incomplete  EOS  for  the  absolute  temperature 

e-e;(v)+ct;(v ) 


T(y,5)  = 


c 


(51) 
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and  pressure 


P( V,  E )  =  -El'  ( V)  -  [E  -  El  (V)]^l , 


the  last  of  which  can  be  rewritten  as 


,  ,  d\nT*(V)  ,  d  E*(V) 

P(V,E)  =  -E - ^-t-T(V) - ff-f 

V  ;  dV  V  JdVT*(V ) 


(52) 


(53) 


4.4  Recovery  of  the  Functions  a(V)  and  MV)  From  the  Quasi-Static  Ambient  Isotherm 
and  the  Hugoniot  Adiabata 

The  potential  S  (V,E)  and  the  kernel  K(V,E )  are  given  by  the  formulas: 


s{v,E)=s-+an^m^E+ 


CT 


V 

C  J  drj- 


2 Eh  (77)  -  IE,  +  P,  [rj  -  V, )  dE *  (77) 


TJ-V, 


dj] 


Eh(tj)~E*t(tj)  +  CT* 


(54) 


and 


K(V,E)  = 


C 


E-E*(V)  +  CT *  ' 


(55) 


Respectively,  the  temperature  and  pressure  are  given  by  the  formulas 


T(V,E) 


1 

K(V,E) 


E-K(v)+cr*  _r  |  g-g;(y) 

c  c 


(56) 


and 


P{V,E) 


dE^{V)  E-E*(V)  +  CT* 

dV  *  Eh(V)-E*  (V)  +  CT* 

E-E*(V)  +  CT*  2EH(V)-2E,  +  P(V -V,) 

v-v.  eh(v)-e*(v)  +  ct * 


(57) 
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4.5  Recovery  of  the  Functions  a(V)  and  A-(V)  From  the  Hugoniot  Energy  and 
Temperature  Measurements 


By  using  equation  31,  we  get  the  potential 

,  .  ,  E-Eh(V)  +  CTh(V ) 

S(V,E)  =  S  +  Cln - ; 


CT„(V) 


V 

\dr! 


+ 


Th(i) 


2 Eh  (tj) - 2E,  +  P*(tj-V*)  dEH  (77) 


TJ-V. 


drj 


(58) 


The  absolute  temperature  and  pressure  are  given  by  the  formulas 


T{V,E) 


1 

K„(V,E) 


e-e„(v)+cth(v) 


and 


P(V,E) 


E  -  E^  (y )  +  crH  (V)  2Fh  (V)  -  2E.  +  P  (V  -  K ) 
CT„{V)  v-v.  ' 

E-EH(v)d\_EH(v)-crH(y}\ 
CTh(V )  dv 


(59) 


(60) 


5.  Conclusion 


We  introduced  a  novel  complete  EOS  —  the  APA  -EOS  —  for  hydrocode,  which  is  an  operator 
on  three  functions  of  one  variable  each:  A(V ) ,  P (£) ,  and  A(V ) .  The  model  is  based  on  the 

simple  phenomenological  assumption  that  the  heat  capacity  Cv  at  constant  volume  is  a  function 

of  the  internal  energy  density  E  only.  Readers  are  reminded  that  the  classical  ideal  case  of  EOS 
is  based  on  the  assumption  that  Cv  is  a  function  of  the  absolute  temperature  T  only,  whereas  the 

Mie-Gruneisen  EOS  is  phenomenologically  equivalent  to  the  assumption  that  Cv  is  a  function  of 
the  entropy  density  S  only  (for  further  information,  refer  to  the  publications  [5-iO]). 

We  established  the  explicit  formulas  permitting  recovery  of  the  function  A(V),  P(£) ,  and 
A(V)  from  different  experimental  measurements. 
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